Electronic and hole spectra of layered systems 
of cylindrical rod arrays: solar cell application 
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We have computed the electronic and hole spectra of a 3D superlattice consisting of layers of 
GaAs rods of finite height arranged in a hexagonal lattice and embedded in an AldGai_dAs matrix, 
alternating with spacer layers of homogeneous AlAs. The spectra are calculated in the envelope 
function approximation, with both light-hole and heavy-hole subbands and hole spin degeneracy 
taken into account. The application of thick spacers allows to investigate the band structure of 
isolated layers of cylindrical rods. We estimate the ultimate efficiency of solar energy conversion in 
a solar cell based on an array of cylindrical quantum dots versus the dot height, and determine the 
optimal value of this parameter. 
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I. INTRODUCTION 

Systems with intermediate energy bands between the 
valence band and the conduction band offer an increased 
efficiency of conversion of solar energy into electricity 1 . 
The theoretical ultimate efficiency limit in solar cells with 
a single intermediate^ band is about 68%^, much above 
the ultimate efficiency limit, estimated at 42%£ in sys- 
tems without intermediate band. 

One of the methods used for generating extra energy 
bands consists in introducing an additional periodicity 
exceeding that of the crystal lattice. This is used in semi- 
conductor superlattices^^—, in which the system is a 
heterostructure with dots, rods or layers of a semiconduc- 
tor material arranged periodically in a matrix of another 
material. If the period of the superlattice exceeds sub- 
stantially that of the crystal lattice, the envelope function 
approximation is applicable. In this case, the parabolic 
conduction band bottom and valence band top are as- 
sumed to be split into a number of energy minibands. 

In ALjGai-dAs, an alloy of zincblende structure, the 
valence band top is composed of bands of light and heavy 
holes, in which the spin degeneracy is lifted for certain 
nonzero wave vectors (k z ^ 0). Thus, the structure of 
minibands within the valence band is rather complesi^. 
Since the values of effective mass of light and heavy holes 
are relatively high compared to those of electrons, mini- 
bands and minigaps in the valence band are much nar- 
rower than the gap between the conduction band and the 
valence band. The structure of minibands in the conduc- 
tion band is more transparent because of the occurrence 
of only one electronic band (in a monolithic material) 
split into relatively large minigaps and minibands in the 
superlattice. 

For typical values of lattice constant of the superlattice 
the lowest conduction miniband, lying below the matrix 
material potential (i.e. in the potential well generated 
by the dots or rods), is well detached from the higher 



minibands, which tend to overlap. In our estimation of 
solar cell efficiency the lowest conduction miniband is as- 
sumed to act as an intermediate band (IB), while the 
minibands in the conduction band (CB) above it, as well 
as the valence band (VB) as a whole, form continuous 
blocks. The ultimate efficiency determined in this study 
refers to this type of band structure. The ultimate ef- 
ficiency reflects how the band structure of the solar cell 
matches the solar spectrum, with inevitable losses in the 
conversion of light quanta into electricity due to the in- 
complete absorption of photons (quanta of energy below 
the width of the narrowest bandgap are not absorbed in 
the system) and thermalization processes (a part of en- 
ergy of photons above the bandgap width is dissipated in 
thermal contact with the lattice). 

The objective of this study is to determine the elec- 
tronic and hole spectra for an isolated array of cylindrical 
quantum dots arranged in a hexagonal lattice. We inves- 
tigate the superlattice spectrum versus the quantum dot 
layer thickness (the height of the cylindrical dots). Next, 
we calculate the ultimate efficiency of a solar cell based 
on the considered structure and examine its dependence 
on the quantum dot layer thickness. 

The paper is organized as follows. After this Introduc- 
tion we discuss the geometry of a 3D superlattice of cylin- 
drical dot arrays. Next, we present the effective mass ap- 
proximation applied to the description of electronic and 
hole states in a heterostructure based on AlGaAs. Our 
results, showing spectra of quantum dot arrays and the 
ultimate efficiency versus the array thickness, are dis- 
cussed in a separate section. This is followed by Conclu- 
sions, which close the paper. 



II. SUPERLATTICE GEOMETRY 

Let us consider a 3D lattice composed of layers of cylin- 
drical quantum dots arranged in a hexagonal lattice, as 
shown in Fig[TJ Acting as quantum wells, the quantum 
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Figure 1. Structure of three-dimensional semiconductor su- 
perlattice in cross section (a) perpendicular and (b) parallel to 
rod axes. Thick arrows in (a) and (b) indicate the position of 
cross section planes. Arrays of rods (bleack area, material A) 
embedded in matrix (gray area, material B) are separated by 
spacer layers (white area, material C). Dashed lines indicate 
unit cell limits (in real space); the first BZ and the irreducible 
part of it, is presented in (c). R is a cylindrical rod radius, ai 
and 3,2 are hexagonal lattice vectors (|ai| = |ai| = a ). 



dots (made of GaAs) are embedded in a matrix material 
(Alo.35Gao.65As), which represents a potential barrier for 
both electrons and holes. The layers of quantum dots 
(quantum dot arrays) alternate with homogenous spacer 
layers of AlAs. Thick spacer layers having a high poten- 
tial barrier of AlAs imply the isolation of individual dot 
layers. 



III. MODEL 

We use the effective mass approximation to calculate 
the electronic and hole states, and assume the electronic 
bands do not interact with the bands of light and heavy 
holes. The calculation of the electronic spectrum is based 
on the Ben Daniel Duke hamiltonian with space- variable 
effective mass to* and variable position of the conduction 
band bottom En^: 
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+ E c (r)^ e (r) =E* e (rll) 

where a = h/2m e , m e denotes the free electron mass and 
"Je is an envelope of the electron function. The determi- 
nation of hole states requires taking into account both 
the light-hole and heavy-hole bands. Since, in general, 
the hole spin degeneracy is lifted in a 3D system, four 
components of the hole envelope function must be 
taken into consideration. Thus, the Schrodinger equa- 



tion for the envelope function become! 
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where 



*h(r) = ( tyfctfr), Vihl(r), * Wli (r), * WlT (r) ) . (3) 

The indexes: hh and Ih denote the components of enve- 
lope function for heavy holes and light holes, respectively. 
The symbols | and J. distinguish bands related to oppo- 
site z-components of spin. The elements: P,Q,R,S of 
matrix in Eq. [2] are given by the following formulae: 
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Luttinger parameters 71, 72, 73 describe the non- 
isotropic effective mass of light and heavy holes, and Ey 
is the position of the valence band top. 

Figure [2] shows the spectrum of the 3D superlat- 
tice calculated by the plane-wave method. For the pe- 
riodic material parameters: to*, Ec, 71, 72,73, Ey, 
Fourier transforms have been calculated analytically. 
The number of reciprocal lattice vectors in the expan- 
sion in a series of plane waves is set so as to obtain 
a satisfactory convergence of results. The dispersion 
relation is determined along the high-symmetry paths 
r-K-M-T-A-H-L-A, K - H and M-L in 
the first Brillouin zone, shown in Fig. [TJ: . 

In the calculation of all the electronic and hole spec- 
tra we have assumed the following values of material pa- 
rameters in the dots (d — 0), the matrix (d — 0.35) 
and the spacers (d = m* = 0.067 + 0.083d, 

E c = 0.944d, 71 = 6.85 - 3.40d, 72 = 2.10 - 1.42d, 
73 = 2.90 - 1.61d, E v = 1.519 - 0.75d. The filling frac- 
tion of dots in the matrix in the quantum dot layer was 
fixed at / = 27ri? 2 /v / 3a 2 = 0.3, and the lattice constant 
of the hexagonal superlattice at a — 50A. The dot layer 
thickness hi and the spacer layer thickness 1h\ were ad- 
justed so as to attain the limiting case of either isolated 
layers (2hi h^) or infinite cylinders (0 ~ 2hi <§C hi). 
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For quantum dots of typical size (see Fig. [2k) and 
layers of typical thickness only the lowest electron mini- 
band is distinctly detached from the other minibands. 
The occurrence of this miniband as an intermediate ex- 
tra level taking part in optical transitions between the va- 
lence band and the block of higher conduction minibands 
increases the efficiency of photon absorption. Besides the 
direct transition, available also in monolithic materials, 
between the valence band and the conduction band (VB- 
CB), a cascade transition valence band - intermediate 
band - conduction band (VB-IB-CB) becomes available 
as well. The power utilized by the photoconverter can be 
evaluated as: 
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where 



I(E 1 ,E 2 ) = 2Tr(k B T s ) 3 /h 3 c 2 f 
£{E) = E/k B T f 



(6) 



denotes the flux of photons of energy ranging from E\ to 
Ei ; E G and Ej are the width of the gap between VB and 
CB and the distance between the IB bottom and the CB 
top, respectively; k B is the Boltzmann constant, and T5 
is the temperature corresponding to the maximum of the 
solar spectrum. The ultimate efficiency of the solar cell 
is defined as the ratio of the utilized power P ou t to the 
power Pi n of the incoming photon flux: 
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where 



P in = 2TT i> (k B T s )yi5h 3 c 2 , 



(7) 



(8) 



denotes the power of the flux of photons emitted by a 
black body at temperature T$- 



IV. RESULTS 

Figure [2] shows the electronic spectra of the quantum 
dot superlattice depicted in Fig. [IJ calculated for differ- 
ent values of thickness /12 of the dot layers and thickness 
2h\ of the spacer layers. The spectra shown in Fig. [2k 
and b correspond to systems with relatively thick spacer 
layers (2h\ — 100A), which implies the dot layers are iso- 
lated from one another. In this case electrons are bound 
within the layer and only allowed to propagate in the x—y 
plane. Along the paths T — A, K — H and M — L (propa- 
gation in the z direction) the relation E(k) is nondisper- 
sive due to the restriction (binding) of the electron mo- 
tion in the layer. Dispersion only occurs when the elec- 
tron energy exceeds the barrier generated by the spacer 
(the barrier level is represented as a horizontal dashed 
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Figure 2. The electronic spectra of the system presented in 
Fig. [T] with GaAs rods (material A) and Alo.35Gao.65 As ma- 
trix (material B), for filling fraction / = 0.3 and lattice con- 
stant a — 50 A in the quantum dot layers. The spectra in (a) 
and (b) correspond to a system with almost isolated layers 
(spacer width 2h\ = 100 A). The dispersive branches corre- 
sponding to the propagation along the rod axis are flat below 
the barriers in spacers (made of AlAs, material C). The higher 
number of modes resulting from a denser quantization along 
the rod axis (for longer rods, (b)) causes minigaps to close. 
The limit of infinite rods (c) is approached by using long rods 
and thin spacers (2hi — 0.2A). The parabolic dispersion re- 
lation is seen to fold in the first BZ for free motion along the 
rod axis. Red dashed line represents the dispersion relation 
in the 2D rod array model (when the motion along z-axis is 
not included. 



line in Fig. [2]) . In a single isolated layer the electron mo- 
tion in the layer plane (the x—y plane) and that in the di- 
rection perpendicular to it (the z direction) are nearly in- 
dependent, as a consequence of the separability of the ef- 
fective potential: E G (x, y, z) = E G>x - y (x, y) + E c , z (z) 



4 



(a strict separation of variables in ([TJ requires the sep- 
aration of variables also in the inverse effective mas*^) . 
Hence, the total electron energy is, approximately, the 
sum of the energy E x - y of the electron motion in the 2D 
periodic potential in the layer and the energy E z of the 
electron motion bound in a ID well of width equal to the 
layer thickness. This is illustrated in Fig. [2j in which 
individual minibands are seen to form through a shift of 
minibands (represented as red dashed lines) correspond- 
ing to the motion in the 2D potential Ec, x -y m the layer 
by the energy levels in the ID potential well Ec,z{z) (cor- 
responding to the position of non-dispersive minibands 
on the path T — A, where k z ^ and k x _ y — 0). The ef- 
fect is illustrated in Fig. 2b, in which the shift of the first 
and the second 2D minibands is indicated by dashed and 
dotted arrows, respectively; solid arrows indicate the en- 
ergy levels E z corresponding to the non-dispersive mini- 
bands on the path V — A. 

The spectra presented in Fig. [2^ and b correspond to 
structures of different quantum dot array thickness. As 
the height of dots increases, the number of bound states 
in the system grows as well. This translates into an in- 
creased number of minibands, and their denser coverage 
of the same energy range. As a result, for dot arrays thick 
enough (Fig. EJd) all the electronic minibands, including 
the lowest one, form a continuous block. For the consid- 
ered structure to provide a basis of a high-efficiency solar 
cell, the dot array thickness should be limited so that the 
minigaps do not close. 

Figure EJ: shows the superlattice spectrum for strongly 
interacting dot arrays. As the spacers are very thin 
(2hi = 0.2A), the system can be regarded, with quite 
good approximation, as a superlattice with cylindrical 
rods of infinite length. In a 2D effective potential gener- 
ated by this type of structure the electron motion along 
the z direction is free. This is evidenced by the parabolic 
dispersion relation along paths T — A, K — H and L — M, 
artificially folded to the first BZ. In an unfolded BZ the 
miniband structure for k z = will be identical with that 
of the 2D superlattice (red dashed line). 

In order to estimate the ultimate efficiency of the solar 
cell we must know the distance Ej between the bottom 
of the intermediate band (the first conduction miniband) 
and the bottom of the conduction band (the continuum 
formed by the minibands above the matrix potential) , as 
well as the gap Eq between the valence band and the 
conduction band (see Fig. [2^). We have determined 
both by investigating the position of absolute minigaps 
and minibands versus dot array thickness hi (see Fig. 
[3^). The calculations were performed for large spacer 
thickness, 2h\ — 100A, ensuring a sufficient isolation of 
the quantum dot arrays. The resulting plot of ultimate 
efficiency versus dot array thickness is shown in Fig. [3Jd. 

All the minibands are easily seen to shift towards 
higher energies with decreasing /12. This is caused by 
an increase in energy component E z related to the quan- 
tization in the z direction. The efficiency of the solar 
cell is in this case mainly limited by the small shift Ej 
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Figure 3. (a) The absolute electronic energy bands (the gray 
area) and gaps (white area) versus rod length hi. (b) The 
ultimate efficiency versus /12, for parameter values adjusted 
so as to apply the isolated layer approximation. The insets in 
(b) present the structures in cross section along rod axis. 

of the intermediate band with respect to the conduction 
band. Thus, changes in width of the gap between the IB 
and the CB translate directly into changes in ultimate 
efficiency of the solar cell. The widest gap corresponds 
to dot array thickness h% ~ 50 A and implies maximum 
ultimate efficiency of the solar cell based on the structure 
under consideration. As evidenced in Fig. [2^, the maxi- 
mum gap width entails the overlap of the two minibands 
originating from the first and the second minibands in 
the 2D dispersion relation. The gap closes, and the IB 
merges into the CB, for h 2 « 72 A, as a result of reduced 
spacing between energy levels E z and the large number 
of minibands in the system. 

V. CONCLUSION 

The thickness of the dot array has a substantial effect 
on the electronic spectrum of the structure considered in 
this study. With dot arrays thick enough, the electronic 
minibands caused by the in-plane periodicity are seen to 
close as a result of a large number of overlapping mini- 
bands. This is due to the dense quantization along the 
direction perpendicular to the layer. The efficiency of a 
solar cell using a quantum dot array, for structures based 
on AlGaAs, is mainly limited by the narrow minigaps in 
the conduction band. The ultimate efficiency decreases 
as the minigaps gradually close as a result of growing 
dot array thickness. The dot array thickness optimal 
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from this point of view corresponds to the widest mini- 
gap; thus, for a hexagonal lattice of GaAs dots embedded 
in Alo.35Gao.65As with filling fraction / = 0.3, the opti- 
mal dot array thickness is h-x — 50A, equal to the lattice 
constant of the 2D heterostructure. 
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